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Abstract. We have determined mass loss rates and gas expansion velocities for a sample of 69 M-type irregular 
(IRV; 22 objects) and semiregular (SRV; 47 objects) AGB-variables using a radiative transfer code to model their 
circumstellar CO radio line emission. We believe that this sample is representative for the mass losing stars of 
this type. The (molecular hydrogen) mass loss rate distribution has a median value of 2.0x lO"'' M© yr~^, and a 
minimum of 2.0x 10~* M0 yr~^ and a maximum of 8xl0~^ M0 yr~^. M-type IRVs and SRVs with a mass loss 
rate in excess of 5xlO~^M0yr^^ must be very rare, and among these mass losing stars the number of sources 
with mass loss rates below a few 10~* Mq yr~^ must be small. We find no significant difference between the IRVs 
and the SRVs in terms of their mass loss characteristics. Among the SRVs the mass loss rate shows no dependence 
on the period. Likewise the mass loss rate shows no correlation with the stellar temperature. 
The gas expansion velocity distribution has a median of 7.0 km s~^ , and a minimum of 2.2 km s^^ and a maximum 
of 14.4 km s~^. No doubt, these objects sample the low gas expansion velocity end of AGB winds. The fraction of 
objects with low gas expansion velocities is very high, about 30% have velocities lower than 5kms~^, and there 
are objects with velocities lower than 3kms~^: V584 Aql, T Ari, BI Car, RX Lac, and L^ Pup. The mass loss rate 
and the gas expansion velocity correlate well, a result in line with theoretical predictions for an optically thin, 
dust-driven wind. 

In general, the model produces line profiles which acceptably fit the observed ones. An exceptional case is R Dor, 
where the high-quality, observed line profiles are essentially flat-topped, while the model ones are sharply double- 
peaked. 

The sample contains four sources with distinctly double-component CO line profiles, i.e., a narrow feature centered 
on a broader feature: EP Aqr, RV Boo, X Her, and SV Psc. We have modelled the two components separately 
for each star and derive mass loss rates and gas expansion velocities. 

We have compared the results of this M-star sample with a similar C-star sample analysed in the same way. The 
mass loss rate characteristics are very similar for the two samples. On the contrary, the gas expansion velocity 
distributions are clearly different. In particular, the number of low-velocity sources is much higher in the M-star 
sample. We found no example of the sharply double-peaked CO line profile, which is evidence of a large, detached 
CO-shell, among the M-stars. About 10% of the C-stars show this phenomenon. 
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1. Introduction 

It has been firmly established that mass loss from the sur- 
face is a very important process during the final stellar 
evolution of low- and intermediate-mass stars, i.e., on the 
asymptotic giant branch (AGB). The mass loss seems to 
occur irrespective of the chemistry (C/0<1 or >1) or the 
variability pattern (irregular, semi-regular, or regular) of 
the star. Beyond these general conclusions the situation 
becomes more uncertain ( lObfe son I999] ). Of importance 
for comparison with mass loss models and for the under- 



standing of AGB-stars in a broader context (e.g., their 
contribution to the chemical evolution of galaxies) is to 
establish the mass loss rate dependence on stellar param- 
eters, such as main sequence mass, luminosity, effective 
temperature, pulsational pattern, metallicity, etc., and its 
evolution with time for individual sources. A crude picture 
has emerged where the average mass loss rate increases 
as the star evolves along the AGB, and where the final 
mass loss rate reached increases with the main sequence 
mass. In addition, there is evidence of time variable mass 
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loss ( Hale et al. 1997 ; Mauron fc Huggins 2000 ), and even 
highly episodic mass loss (Olofsson et al. 2000). There 
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is a dependence on luminosity in the expected way, i.e., 
an increase with increasing luminosity, but it is uncertain 
how strong it is. The same applies to the effective tem- 
perature where a decrease with increasing temperature is 
expected. Regular pulsators clearly have higher mass loss 
rates than stars with less regular pulsation patterns. The 
dependence of the total mass loss rate on metallicity ap- 
pears to be weak, but the dust mass loss rate decreases 



with decreasing metallicity. See Habing (1996) for a sum- 
mary of evidence in favour of this general outline. 

The mechanism behind the mass loss remains un- 
known, even though there are strong arguments in favour 
of a wind which is basically pulsation-driven, and where 
the highest mass loss rates and gas expansion velocities 
are reached through the addition of radiation pressure on 
dust (iHofner fc Dorfi 1997j ; [Winters et al. 2000] ). A way 
to study this problem is to use samples of low mass loss 
rate stars for which stellar parameters can be reasonably 
estimated using traditional methods. These samples also 
contain objects with quite varying pulsational character- 
istics, and has, as it turned out, quite varying circumstel- 
lar characteristics. Olofsson et al. (1993| ) presented such 
a study of low mass loss rate C-stars using CO multi- 
transition radio data. These data were subsequently anal- 



ysed ir more detail by Schoier & Olofsson (2001) using a 



radiative transfer model. In the same spirit Kcrschbaum 
& Olofsson (1999) presented a major survey of CO radio 
line emission from irregularly variable (IRV) and semireg- 
ular ly variable (SRV) M-type AGB-stars. They increased 
the number of IRVs (22 detections), in particular, and 
SRVs (43 detections) detected in circumstellar CO emis- 
sion substantially («60% of the SRVs and all but one of 
the IRVs were detected for the first time). Young (1995) 
and Groenewegen et al. (1999) have made extensive sur- 
veys of short-period M-Miras. 

In this paper we use the radiative transfer method of 



Schoier fc Olofsson (2001[) to estimate reasonably accu- 



rate mass loss rates and gas e xpansion velocities for the 
Kerschbaum fc Olofsson (1999 ) sample. Comparisons be- 



tween these properties and other stellar properties are 
done, as well as comparisons with the results for the C-star 
sample. 

2. The sample 

The original sample consisted of all IRV and SRV AGB- 
stars of spectral type M (determined by spectral classifi- 
cation, or using the IRAS LRS spectr a) in the Genera l 
Catalogue of Variable Stars [GOVS; [Kholopov (l"990| )] 
with an IRAS quality flag 3 in the 12, 25, and 60 /im bands. 
From this sample we selected for the CO radio line ob- 
servations objects with IRAS 60/im fluxes, 5*60, typically 
>3 Jy. Since Sqq for a star with a luminosity of 4000 L© 
(see below) and a temperature of 2500 K is 34 Jy, 6Jy, 
and 1.4 Jy at a distance of 100 pc, 250 pc, and 500 pc, re- 
spectively, we added a colour selection criterium. Only 
stars redder than 1.2 mag in the IRAS [12]-[25] colour 
were observed, thus biasing the sample towards stars with 



detectable circumstellar dust envelopes. There is a possi- 
bility that stars with detectable gas mass loss rates, but 
with very little circumstellar dust, were missed due to this. 
About 50% of the objects, i.e., 109 sources, were subse- 
quently searched for circumstellar radio line emission. 

The CO (J=1^0, 2^1, 3^2, and 4^3) data which 
are used as the observational constraints for the mass 
loss rate determinations in this paper were presented in 



Kerschbaum fc Olofsson (1999). They were obtained using 



the 20 m telescope at Onsala Space Observatory (OSO), 
Sweden, the 15 m Swedish-ESO Submilhmetre Telescope 
(SEST) on La Silla, Chile, the IRAM 30 m telescope 
on Pico Veleta, Spain, and the James Clerk Maxwell 
Telescope on Mauna Kea, Hawaii. A few additional sources 
were observed at OSO in May 2000, see Sect. ||. In to- 
tal, 69 stars were detected, 22 IRVs and 47 SRVs, and 
45 were detected in at least two transitions, and 20 in 
at least three transitions. The detection rate was rather 
high, about 60%, and the conclusions drawn in this paper 
should be representative for the stars in our sample. 

The distances, presented in Table |[ were estimated 
using an assumed bolometric luminosity of 4000 Lq . This 
value was chosen in agreement with the typical values 



derived by Kerschbaum et al. (1997) and Barthes et al 
(1999| ) for objects with similar properties. Morover, the 



13 objects in our sample having Hipparcos parallax errors 
better than 20% have a mean luminosity of 4200 Lq (with 
a standard deviation of 1900 Lq, i.e., consistent with the 
assumed luminosity and the parallax uncertainty). For a 
statistical study of a large sample of stars these distance 
estimates are adequate (and they were used also for stars 
with reliable Hipparcos distances to avoid systematic dif- 
ferences) , although the distance estimate for an individual 
star has a rather large uncertainty. 

The apparent bolometric fluxes were obtained by in- 
tegrating the spectral energy distributions ranging from 
t he visual data over the near-infrared to the I RAS-range 
(Kerschbaum fc Hron 1996; Kerschbaum 199£). 



3. Observations 

A few additional sources were observed at OSO in 
May 2000 with the same instrumental setup as used by 



Kerschbaum fc Olofsson (1999 ). Relevant information on 



the instrumental setup, and the method used to derive the 
line proflle properties and the upper limits can be found 
in Kcrschbaum fc Olofsson (1999 ). A summary of the ob- 
servational results are given in Table (where we give 
the velocity- integrated intensities, /mb, and antenna tem- 
peratures, Tmb, in main beam brightness scale, the stellar 
velocities as heliocentric, t^hoi, and Local Standard of Rest, 
't'LSR, velocities, and the gas expansion velocity, Vc), and 
the detections are presented in Fig. |^. 

4. Radiative transfer 

We have chosen to use the Monte Carlo method to de- 
termine the excitation of the CO molecules as a function 



Fig. 1. CO(J=1->0) detections at OSO in May 2000 



of distance from the central star ( Bernes 1979| ). This is 
a versatile method, close to the physics, which has been 
applied to CO radio and far-IR line emission (prosas fc 
Mente]|i 1997|; |Ryde et al. 1999]; ^choier ct al. 2002|), CO 



near-IR line scattering ( Ryde fc Schoier 2001|), and to the 



OH 1612MHz maser emission (van Langevelde & Spaans 
1993|3Sch6ier & Olofsson (2000, 2001) apphed it to CO 
radio line data of samples of AGB carbon stars, and the 
results suggest that the modelling of circumstellar CO 
emission is one of the most reliable methods for estimat- 
ing mass loss rates. The code has been benchmark tested 
together with a number of other radiative transfer codes 
( van Zadelhoff ct al. 2002 ). Below we give a brief intro- 
duction to the modelling work, while for details we refer 
to 



Schoier fc Olofsson (2001) 



4.1. The CO molecule 

In the excitation analysis of the CO molecule we used 
40 rotational levels in each of the ground and first ex- 
cited states. Excitation to the v=2 state can be ignored, 
since the v=l state is not well populated. The radia- 
tive transition probabilities and energy levels were taken 
from [Chandra et al. (1996| ). The coUisional rate coeffi- 



cients (CO— para-H2) for rotational transitions are based 
on the results in Flower fc Launay (1985| ). These are fur- 
ther extrapolated for J>11 and for temperatures higher 
than 250 K. We neglect collisional transitions between the 
vibrational states because of the low densities and the rel- 
atively low temperatures. 



Recently, Flower (2001) presented revised and ex- 
tended collisional rates for CO— H2. Individual rates are 
generally different from those previously published in 



Flower fc Launay (1985), with discrepancies as large as 



a factor of two in some cases. In addition, for temper- 
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atures above 400 K the rates from Schinke et al. (1985) 
were used and further extrapolated to include transitions 
up to J=40. To test the effects of the adopted set of col- 
lisional rates a number of test cases with the new rates 
assuming an H2 ortho-to-para ratio of three were run. We 
found that for the relatively low mass loss rate stars of 



interes ; here, where excitation of CO from radiation dom- 



inates over that from collisions with H2, the adopted set 
of coUisional rates is only of minor importance. 

4.2. The circumstellar model 

4.2.1. The geometry and kinematics 

We have adopted a relatively simple, yet realistic, model 
for the geometry and kinematics of the CSE. It is as- 
sumed that the mass loss is isotropic and constant with 
time. As will be shown below the median mass loss rate 
is 2xlO~^Moyr~^ and the median gas expansion veloc- 
ity is 7.0 km s^^ for our sample. This results in an extent 
of the CO envelope of about 2x10"^^ cm (see Sect. 4.2. 3] ), 
and the corresponding time scale is about 1000 yr. There 
is now some evidence for mass loss rate modulations of 



the thermal width of CO is about 0.3(Tk/100)"-5 kms"!]. 
This can be an important parameter for low mass loss 
rate objects since it affects the radial optical depths and 
hence the effectiveness of the radiative excitation. The 
constraints on this parameter are rather poor. The most 
thorough analysis in this connection is the one by Huggins 
Healy (1986| ) . They modelled in detail the circumstellar 



AGB-stars on this tim e scale ( Mauron fc Huggins 200C ; 
Marengo et al. 2001 ), and therefore the assumption of 
a constant mass loss rate may be questionable, and this 
must be kept in mind when interpreting the results. The 
gas expansion velocity, Vc, is assumed to be constant with 
radius. This ignores the more complex situation in the in- 
ner part of the CSE, but the emission in the CO lines 
detected by radio telescopes mainly comes from the ex- 
ternal part. We further assume that the hydrogen is in 
molecular form in the region probed by the CO emission 
( |Glassgold fc Huggins 19*83 ). As a consequence of these 
assumptions the H2 number density follows an r^^-law. 

In the case of low mass loss rate objects the inner ra- 
dius of the CSE will have an effect on the model intensi- 
ties. The reason is that radiative excitation plays a role in 
this case and the absorption of pump photons at 4.6 /im 
depends (sensitively) on this choice. Wc set the inner ra- 
dius to Ix lO^"* cm (~3 R*), i.e., generally beyond both the 
sonic point and the dust condensation radius. The uncer- 
tainty in the mass loss rate estimate introduced by this 
assumption is discussed in Sect. O, Strictly, speaking the 
assumption of a constant expansion velocity from this in- 
ner radius is very likely not correct. An acceleration region 
will enhance the radiative excitation and hence may have 
an effect on the estimated mass loss rate. However, as will 



be shown in Sect. 4.2, the dependence on the inner ra- 
dius is rather modest, suggesting that the properties of 
the inner CSE are not crucial for the mass loss rate de- 
termination. We have therefore refrained from introducing 
yet another parameter, i.e., a velocity law parameter. 

In addition to thermal broadening of the lines micro- 
turbulent motions contribute to the Doppler broadening of 
the local line width. We assume a turbulent velocity width, 
Vt, of 0.5kms^^ throughout the entire CSE, i.e., the same 



CO line self-absorption in the high mass loss rate car- 
bon star IRC-f 10216 and derived a value of 0.9 km s^^. In 



Sect. 4.3 we discuss the uncertainty in the mass loss rate 
estimate introduced by this parameter. 

4.2.2. Heating and cooling processes 

We determine the kinetic gas temperature in the CSE by 
taking into account a num ber of heating and cooling pro- 
cesses ( Groencwcgcn 1994 ). The primary heating process 
is the viscous heating due to the dust streaming through 
the gas medium. A drift velocity between the gas and 
the dust is calculated assuming a dust-driven wind (Kwok 
1975 ), but for the low mass loss rate stars in this study 
the radiation pressure on dust may not be very efficient, 
i.e., the driving of the gas may be due to something else. 
However, as will be explained below, the gas-dust heating 
term is nevertheless very uncertain, and we use it as a 
free parameter in our model. Additional heating is due to 
the photoelectric effect, i.e., he ating by electrons ej ected 
from the grains by cosmic rays (Huggins et al. 1988), but 
for our low mass loss rate stars this has a negligible effect 
inside the CO envelope. 

There are three major cooling processes, adiabatic ex- 
pansion of the gas, CO line cooling, and H2O hue coohng. 
The CO line cooling is taken care of self-consistently by 
calculating it s magnitude after each iteration using the 
expression of Crosas fc Menten (1997). H2O line cooling 



is estimated using the results from Neufeld fc Kaufman 



(1993 ). They calculated the H2O excitation using an es- 



cape probability method and estimated the radiative cool- 
ing rates for a wide range of densities and temperatures. 
The H2O abundance is set to 2x10"'* and the envelope 
sizes used are based on the results of Netzcr fc Knapp 



(1987). In addition, H2 line cooling is taken into account 



(Groencwegcn 1994), but this has negligible effect in the 
regions of interest here. 

When solving the energy balance equation a number 
of (uncertain) pa rameters describing the du st are intro- 
duced. Following Schoicr fc Olofsson (2001 ) we assume 
that the Qp.F-para-meter, i.e., the flux-averaged momen- 
tum transfe r effeciency from th e dust to the gas, is equal 
to 0.03 [see Habing et al. (1994 ) for details], and define a 
new parameter which contains the other dust parameters. 



0.01 



2.0 gcm- 



0.05 ^m' 



(1) 



value as used by Schoier fc Olofsson (2001) [for reference 



where ip is the dust-to-gas mass ratio, pgi- the dust grain 
density, and Ogr its radius. The normalized values are the 
ones used to fit the CO radio line emission of IRC-l-10216 
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using this model ( ^choier fc Olofsson 2001 ), i.e., h=l for 
this object. 

4.2.3. The CO fractional abundance distribution 

We assume that the initial fractional abundance of CO 
with respect to H2, fco, is 2xl0~"^, w hich is the same 
value as used by Kahanc fc Jura (1994 ) in their analysis 
of CO radio line emission from M-stars. This is essentially 
a free parameter, although its upper limit is given by the 
abundance of C (i.e., 7x10"'' is the upper limit for a solar 
C abundance). Due to photodissociation by the interstellar 
radiation field the CO abundance starts to decline rapidly 
at a radius, which, for not too low mass loss rates, depends 
on the mass loss rate. Calculations, taking into account 
dust-, self- and H2-shielding, and chemical fractionation, 
have been performed by Mamon et al. (1988|) and Doty 
fc Leung (199^ ). Here we use the results of Mamon et al 
(1988| ) in theway adopted by [Schoier fc Olofsson (200l| )". 

An approximate expression for the photodissociation 
radius, Tp, consists of two terms, the unshielded size due to 
the expansion, which is independent of the mass loss rate, 
and the size due to the self-shielding, which scales roughly 
as (/coM)°-5 dSchoier fc Olofsson 200l| ). These terms are 
equal at a mass loss rate of about 4x 10~^('yc/7)^ M0 yr^^ 
for the adopted CO abundance {vq is given in kms^^). 
That is, self-shielding plays a role for essentially all of our 
objects. We note that for low mass loss rate objects the 
spatial extent of the CO envelope is particularly important 
since the spatial extent of the CO line emission is limited 
by this, and not by excitation. Sect. 4.3. 



4.2.4. The radiation fields 

The radiation field is provided by two sources. The cen- 
tral radiation emanates from the star, and was estimated 
from a fit to the spectral energy distribution (SED), usu- 
ally by assuming two blackbodies, one representing the di- 
rect stellar radiation and one the dust-processed radiation 



(Kerschbaum fc Hron 1996). The dust mass loss rates of 



our sample stars are low enough that the latter can be ig- 
nored. The temperatures derived are given in Table ^. The 
stellar blackbody temperature Tbb derived in this manner 
is generally about 500 K lower than the effective temper- 



ature of the star (Kerschbaum fc Hron 1996). The second 



radiation field is provided by the cosmic microwave back- 
ground radiation at 2.7 K. 

4.3. Dependence on parameters 

We have checked the sensitivity of the calculated in- 
tensities on the assumed parameters for a set of model 
stars. The model stars have nominal mass loss rate 
and gas expansion velocity combinations which are 
characteristic of our sample: (4x 10~* yr^^, 5kms^^), 
(IxlQ-^Moyr-i, Ikms'^), and (SxIQ-^Mq yr"!, 
lOkms"^). They are placed at a distance of 250pc (a 



typical distance of our stars), and the nominal values 
of the other parameters are L=4000Lq, Tbb=2500K, 
h=0.2 (the value adopted for the majority of our stars, 
Sect. ^?^ ), ri=2xl0^^cm (this is twice the inner radius 
used in the modelhng), t!t=0.5kms^^, /co=2xl0^^, 
and Tp is calculated from the photodissociation model 



(Sect. |T|). The CO(J=1^0), C0(J=2^1), and 
CO(J=3^2) lines are observed with beam widths of 
33", 23", and 14", respectively. These are characteristic 
angular resolutions of our observations. Note that, to 
some extent, the presented results are dependent on 
the assumed angular resolution since resolution effects 
may play a role. We change all parameters (except the 
expansion velocity) by —50% and +100% and calculate 
the velocity-integrated intensities. 

The results are summarized in Table ^ in terms of per- 
centage changes. Although the dependences are somewhat 
complicated there are some general trends. The line in- 
tensities are sensitive to changes in the mass loss rate, 
more the lower the mass loss rate, and hence are sensi- 
tive measures of this property. There is a dependence on 
luminosity, in particular for low- J lines for low mass loss 
rates and for high- J lines for higher mass loss rates. The 
dependence on the (uncertain) /i-parameter is fortunately 
rather weak. The dependence on the photodissociation ra- 
dius is substantial, in particular for the low- J lines and for 
low mass loss rates. The dependence on the inner radius is 
weak, and so is the dependence on the turbulent velocity 
width. Thus, we conclude that for our objects the CO ra- 
dio line intensities are good measures of the mass loss rate, 
but it shall be kept in mind that they are rather dependent 
on the uncertain photodissociation radius and, to some ex- 
tent, on the assumed luminosity. A similar sensitivity anal- 
ysis for C-stars were done by ^choier fc Olofsson (2001 ), 
and studies of the parameter dependence w ere done by 



Kwan fc Webster (1993| ) 



Kastner (19921) and 

There is also a dependence on the adopted CO 
abundance. For the low mass loss rates considered here 
(<5x 10~^ Mq yr"-'^) a constant product fcoM produces 
the same model line intensities. Hence, mass loss rates for 
a different value of /co are easily obtained. The reason 
for this behaviour is that the size of the emitting region is 
photodissociation limited rather than excitation limited, 
i.e., for our objects it scales roughly as {fcoM)^-^, see 



Sect. 4.2.5 . Furthermore, the energy levels are to a large 
extent radiatively excited, i.e., the density plays less of 
a role for the excitation. Therefore, since the emission is 
optically thin, a change in the abundance must be com- 
pensated by an equal, but opposite, change in the mass 
loss rate to keep the calculated intensities unchanged. 

An additional uncertainty is due to the somewhat 
crude treatment of H2O line cooling. The modelling shows 
that this cooling process has an effect on the temperature 
structure of the CSE. However, it is found to be of impor- 
tance only in the inner warm and dense part of the CSE 
where H2O is abundant. The size of the H2O envelope 
is only about a tenth of the CO envelope. This small re- 
gion contributes only marginally to the line intensities of 



Fig. 2. SW Vir C0( J=l->0 and 2-^1) spectra obtained with the SEST, and C0( J=3^2 and 4^3) spectra obtained 
with the JCMT (histograms). The hne profiles from the best-fit model are shown as solid, thin lines (the beam size is 
given in each panel) 



the observed low- J transitions (see for instance the weak 
dependence on ri). 

5. Model results 

5.1. Best-fit model strategy 

The radiative transfer analysis produces model brightness 
distributions. These are convolved with the appropriate 
beams to allow a direct comparison with the observed 
velocity-integrated line intensities and to search for the 
best fit model. There are two remaining parameters to 
vary in this fitting procedure, the mass loss rate M and 
the /i-parameter. These two parameters were varied un- 
til the best-fit model was found. The quality of a fit was 
quantified using the chi-square statistic, 

2 _ 1 [-^mod.i ~ -^obs.i)]^ ,r.\ 

^-d - ^ 2. -2 , (2) 

^ i=l » 



where / is the total integrated line intensity, ct; the uncer- 
tainty in observation i, p the number of free parameters, 
and the summation is done over all independent obser- 
vations N. The errors in the observed intensities are al- 
ways larger than the calibration uncertainty of ^^20%. We 
have chosen to adopt ai = 0.2/obs,i to put equal weight 
on all lines, irrespective of the S/N-ratio. Initially a grid, 
centered on M=10~^ yr~-^ and h=0.1, with step sizes 
of 50% in M and 100% in h was used to locate the x^- 
minimum. The final parameters were obtained by decreas- 
ing the step size to 25% in M and 50% in h and by in- 
terpolating between the grid points. The final chi-square 
values for stars observed in more than one transition are 
given in Table ^. The line profiles were not used to dis- 
criminate between models, but differences between model 
and observed line profiles are discussed in Sect. 

In general, the model results fit rather well the ob- 
served data as can be seen from the chi-square values. For 



Fig. 3. The /i-parameters derived from the radiative 
transfer analysis plotted against the density measure 
M /vc- Objects with three or more transitions observed 
are marked with filled circles, while those with only two 
transitions observed are marked with open circles 



uncertainties are such that we can only conclude that both 
the M-type and C-type CSEs due to low mass loss rates 
appear to have dust properties significantly different from 
those of C-type CSEs due to high mass loss rates (how- 
ever, all mass loss rate determinations are made using the 
same value for the flux-averaged dust momentum trans- 
fer efficiency, which determines the gas-dust drift velcoity 
and hence affects the heating of the CSE, while in reality 
it may depend on the mass loss rate). We can also con- 
clude that the gas-CSEs due to low mass loss rates are 
cooler than expected from a simple extrapolation of the 



results for IRC+10216 (Groenewegen et al. 1998; Crosas 



Menten 1997; 3choier fc Olofsson 2000) 



5.3. The mass loss rates 

The estimated mass loss rates are given in Table |^ 
(rounded off to the number nearest to 1.0, 1.3, 1.5, 2.0, 
2.5, 3, 4, 5, 6, or 8, i.e., these values are separated by 
about 25%), and their distributions are shown in Fig. |j 
(four sources with clearly peculiar line profiles are dis- 



cussed separately in Sect. 6/7). We estimate that, within 



the adopted circumstellar model, the derived mass loss 
rates are uncertain by about ±20% for those sources with 
three or more observed transitions, since the line intensi- 
ties are very sensitive to changes in this parameter (see 
Table The uncertainty increases when less than three 
transitions are observed, generally ±50%, but it may be 
as bad as a factor of a few for objects with low h-values. 
To this should be added the uncertainty due to the dis- 
tance, the luminosity, the photodissociation radius, the 
fractional CO abundance, the coUisional cross sections, 
and the pointing/calibration. Nevertheless, as a whole, we 
believe that these are the most accurate mass loss rates 
determined for these types of objects, but on an absolute 
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Table 3. Source parameters and model results 



Source 


Var. 


P 


Di 


Tbb 


M 






h 




N 




type 


fdaysl 


[pc] 


[Kl 

L J 


[10"'^ M© yr-^l 


[km s"M 

L J 


[10'« cml 








BC And 


Lb 




450 


2510 


2.0 


4.0 


2.0 






1 


CE And 


Lb 




740 


2720 


5 


10.5 


2.5 






1 


RS And 


SRa 


136 


290 


2620 


1.5 


4.4 


1.6 




0.7 


2 


UX And 


SRb 


400 


280 


2240 


4 


12.8 


2.1 




1.9 


2 


TZ Aql 


Lb 




470 


2460 


1.0 


4.8 


1.3 






1 


V584 Aql 


Lb 




390 


2340 


0.5 


2.2 


1.2 






1 


AB Aqr 


Lb 




460 


2580 


1.3 


4.2 


1.5 






1 


SV Aqr 


Lb 




470 


2180 


3 


8.0 


2.1 




9.1 


2 


e Aps 


SRb 


119 


110 


2620 


0.4 


4.5 


0.8 




9.8 


2 


T Ari 


SRa 


317 


310 


2310 


0.4 


2.4 


0.9 






1 


RX Boo 


SRb 


340 


110 


2220 


5 


9.3 


2.6 




1.4 


2 


RV Cam 


SRb 


101 


350 


2570 


2.5 


5.8 


2.0 




0.4 


2 


BI Car 


Lb 




430 


2420 


0.3 


2.2 


0.9 






1 


V744 Cen 


SRb 


90 


200 


2750 


1.3 


5.3 


1.5 


0.05 


6.0 


3 


SS Cep 


SRb 


90 


340 


2580 


6 


10.0 


2.7 




3.3 


2 


UY Get 


SRb 


440 


300 


2400 


2.5 


6.0 


2.0 


0.2 


0.4 


3 


CW Cnc 


Lb 




280 


2400 


5 


8.5 


2.5 


0.25 


3.2 


3 


RY CrB 


SRb 




550 


2340 


4 


5.7 


2.5 






1 


R Crt 


SRb 


160 


170 


2130 


8 


10.6 


3.0 


0.3 


0.7 


4 


AF CvK 


SRb 




300 


2840 


0.8 


3.5 


1.2 






1 


W CvK 


SRb 


131 


130 


2670 


1.0 


8.3 


1.3 




0.7 


2 


U Del 


SRb 


110 


210 


2720 


1.5 


7.5 


1.5 




10.5 


2 


R Dor 


SRb 


338 


45 


2090 


1.3 


6.0 


1.4 


0.7 


1.6 


3 


AH Dra 


SRb 


158 


340 


2680 


0.8 


6.4 


1.1 






1 


CS Dra 


Lb 




370 


2580 


6 


11.6 


2.7 


0.05 


3.2 


3 


S Dra 


SRb 


136 


270 


2230 


4 


9.6 


2.2 


0.3 


0.5 


3 


SZ Dra 


Lb 




510 


2580 


6 


9.6 


2.7 






1 


TY Dra 


Lb 




430 


2300 


6 


9.0 


2.8 




1.0 


2 


UU Dra 


SRb 


120 


320 


2260 


5 


8.0 


2.5 




5.0 


2 


g Her 


SRb 


89 


100 


2700 


1.0 


8.4 


1.3 






1 


AK Hya 


SRb 


75 


210 


2430 


1.0 


4.8 


1.3 


0.15 


1.5 


4 


EY Hya 


SRa 


183 


300 


2400 


2.5 


11.0 


1.8 




1.1 


2 


FK Hya 


Lb 




310 


2630 


0.6 


8.7 


1.0 






1 


FZ Hya 


Lb 




330 


2460 


2.0 


7.8 


1.6 




0.0 


2 


W Hya 


SRa 


361 


65 


2090 


0.7 


6.5 


1.0 






1 


RX Lac 


SRb 




250 


2450 


0.8 


2.2 


1.6 






1 


RW Lcp 


SRa 


150 


400 


2150 


0.5 


4.4 


0.9 






1 


RX Lep 


SRb 


60 


150 


2660 


0.5 


3.5 


1.0 




0.2 


2 


SY Lyr 


SRb 




640 


2410 


6 


4.6 


3.6 






1 


TU Lyr 


Lb 




420 


2470 


3 


7.4 


2.1 




1.3 


2 


U Men 


SRa 


407 


320 


2160 


2.0 


7.2 


1.7 






1 


T Mic 


SRb 


347 


130 


2430 


0.8 


4.8 


1.2 




1.6 


2 


EX Ori 


Lb 




470 


2490 


0.8 


4.2 


1.3 




2.2 


2 


V352 Ori 


Lb 




250 


2560 


0.5 


8.4 


0.9 






1 


S Pav 


SRa 


381 


150 


2190 


0.8 


9.0 


1.1 






1 


SV Peg 


SRb 


145 


190 


2330 


3 


7.5 


2.1 


0.1 


4.2 


3 


TW Pee 


SRb 


929 


200 


2690 


2.5 


9.5 


1.8 


0.25 


3.5 


4 


V PsA 


SRb 


148 


220 


2360 


3 


14.4 


1.9 






1 


I? Pup 


SRb 


141 


85 


2690 


0.2 


2.3 


0.7 


0.05 


0.9 


4 


OT Pup 


Lb 




500 


2630 


5 


9.0 


2.6 




4.8 


2 


Y Scl 


SRb 




330 


2620 


1.3 


5.2 


1.5 




2.0 


2 


CZ Ser 


Lb 




440 


2150 


8 


9.5 


3.2 




0.2 


2 


Ser 


SRb 


100 


170 


2500 


1.5 


14.4 


1.5 






1 


SU Sgr 


SRb 


60 


240 


2090 


4 


9.5 


2.3 




4.7 


2 


UX Sgr 


SRb 


100 


310 


2520 


1.5 


9.5 


1.5 






1 


V1943 SgT 


Lb 




150 


2250 


1.3 


5.4 


1.4 


0.5 


9.2 


3 


V T(.-l 


SRb 


125 


290 


2200 


2.0 


0.8 


i.G 




9.5 


2 



^ Distance derived assuming a luminosity of 4000 Lq 



Fig. 5. Gas expansion velocity distribution for the whole sample (left panel), as well as for the IRVs (middle) and 
SRVs (right) separately. The objects with double-component line profiles are excluded 



scale they are uncertain by at least a factor of a few for an 
individual object. A comparison with the mass loss rates 
estimated by Knapp et al. (19981) for the six stars in com- 



mon shows differences by less than a factor of two, except 
in the case of RT Vir for which we derive a five times 
higher value. Note that the mass loss rates given are not 



corrected for the He-abundance, i.e., they are molecular 
hydrogen mass loss rates. 

The detailed radiative transfer presented here results 
in considerably higher mass loss rates than those obtained 



with the simpler analysis in Kerschbaum et al. (1996) and 



Kerschbaum & Olofsson (1998). For the 43 objects in com- 
mon we derive mass loss rates which are on average ten 



10 
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times higher for the same distances and CO abundance 
(the median difference is six). This confirms the conclu- 
sion by 3choier & Olofsson (2001) that the formulae of 
Knapp fc Morris (1985D and [Kastner (1992D lead to sub- 
stantially underestimated mass loss rates for low mass loss 
rate objects. This is no surprise since both formulae were 
calibrated against IRC-l-10216 (for which h=l), and in the 
case of Knapp fc Morris (1985| ) an older CO photodisso- 
ciation model was used. 



5.4. The gas expansion velocities 

The gas expansion velocities given in Table || are obtained 
in the model fits. Hence, they should be somewhat more 
accurate than the pure line profile fit results given in 
Kerschbaum & Olofsson (1999), since for instance the ef- 
fect of turbulent broadening is taken into account. The 
former are in general somewhat lower than the latter. We 
estimate the uncertainty to be of the order ±10%. The 
uncertainty is dominated by the S/N-ratio since the spec- 
tral resolution is in most cases more than adequate. A 
significant fraction of the sources has gas expansion ve- 
locities lower than 5kms~^, and for these the assumption 
of a turbulent velocity width of 0.5kms~^ will have some 
effect on the expansion velocity estimate. The gas expan- 



lection effect. All our stars are included in the GCVS, and 
even though this favours less obscured stars the GCVS 
contains many stars with mass loss rates in excess of 
this value. Thus, there appears to be an upper limit for 
the mass loss rate of an M-type IRV or SRV on the 
AGB. Regular pulsators of M-type, Miras and OH/IR- 
stars, clearly reach significantly higher mass loss rates, 
and hence the regularity of the pulsation and its ampli- 
tude play an important role for the magnitude of the mass 
loss rate. Some caution must be exercised here though. We 
are averaging the mass loss rate over a time scale of about 
one thousand years, and there are indications that the 
mass loss rates of IRV/SRVs are more time- variable, on 
shorter time scales than this, than are the mass loss rates 
of the Miras ( Marengo et al. 2001 ). This would lead to an, 
on the average, lowered mass loss rate of an IRV/SRV. 

The decrease in the number of objects with low mass 
loss rates could indicate an effect of limited observational 
sensitivity. However, a plot of the mass loss rate as a func- 
tion of the distance suggests that this is not the case. 
Fig. H. We detect low mass loss rate objects out to about 
500 pc, and beyond this only a few higher mass loss rate 
objects are detected. That is, nearby stars with mass 
loss rates below lO~^M0yr~^ should be detectable. [We 



sion ve |locity distribution tor the whole sample, as well as^ 
those ot the IRVs and HRVs separately, are show n m I' lg. p|" 
(excluding the double-component sources, Sect. 6.7). 



checked all the non-detections reported by Kerschbaum & 



Olofsson (1999), and found that in no case do they provide 



6. Discussion 

In this section we present and discuss a number of results 
based on the derived mass loss rates and gas expansion 
velocities. Extensive comparisons are made with the re- 
sul ts of the C-type IRVs and SRVs (39 objects) analysed 
by Schoier fc Olofsson (2001 ) using the same methods as 
in this study. 

6.1. The mass loss rate distribution 

The distribution of the derived mass loss rates have a 
median value of 2.0x 10^^ Mq yr"-'^, and a minimum of 
2.0x10"® M0yr-i and a maximum of 8x 10""^ Mq yr'^ 
There is no significant difference between the IRVs and the 
SRVs, but the number of IRVs is quite low. We believe that 
these mass loss rate distributions are representative for the 
mass losing M-type IRVs and SRVs on the AGB (see argu- 
ments below) . We find no significant difference when com- 
paring with the sample of C-type IRVs and SRVs, where 
the median was 1.6x 10"^ M© yr~^. Therefore, the mass 
loss rates of these types of variables appear independent 
of chemistry (also for the C-stars there is no dependence 
on the C/O-ratio). This conclusion rests on the uncertain 
assumptions of the CO fractional abundance (10~^ in the 
case of the C-star sample, as opposed to the value 2x 10"'' 
used here for the M-stars). 

Most notable is the sharp cut-off at a mass loss rate 
slightly below 10~^ Mq yr~^. This is very likely not a se- 



an upper limit in the mass loss rate which is significantly 
lower than a few 10~® M© yr"-'^.] Hence, we interpret the 
trailing off at low mass loss rates as due to a lack of such 
sources among the mass losing M-type IRVs and SRVs on 
the AGB. However, our sample is limited by the IRAS 
colour [12]-[25] (Sect. ||). Therefore, it is possible that 
there exists M-type AGB-stars with mass loss rates lower 



than our limit of about 10 M 



The case for the C- 



stars is different, pchoier fc Olofsson (2001]) also derived a 
lower limit of about 10~® M© yr~^, but this is based on a 
K-magnitude limited sample, for which the K-magnitude 
is expected to be relatively constant, where all stars within 
about 500 pc were detected ( Olofsson et al. 1993 ). 

We have also compared the derived mass loss rates 
with the periods of the SRVs, Fig. |^. The conclusion by 
Kerschbaum et al. (1996) that for these objects the period 
of pulsation plays no role for the mass loss rate still holds. 
The apparent absence of stars with periods in the range 
200-300 days is most probably due to a distinct division 
into two pulsational modes, one operating only below 200 
days and one only above 300 days. Likewise, for the C- 
SRVs we find no correlation between mass loss rate and 
period (the gap between 200 and 300 days does not exist 
for these stars). 

6.2. Mass loss and stellar temperature 

It has turned out to be very difficult to derive the mass 
loss rate of an AGB-star from first principles. Some at- 
tempts have nevertheless been made and they all indicate 
a strong dependence on the stellar effective temperature, 



Fig. 8. The derived mass loss rate as a function of the 
stellar blackbody temperature 



6.3. The gas expansion velocity distribution 

The gas expansion velocities derived from the model fits 
have a distribution with a median for the whole sample 
of 7.0 km s^^, and a minimum of 2.2 km s^^ and a maxi- 
mum of 14.4 km s~^, i.e., clearly these objects sample the 
low gas expansion velocity end of AGB winds. Similar re- 



sults have been obtained for sho rt-period M-Miras ( Young 
1995 ; Groenewegen et al. 1999| ). We find no apparent dif- 



ference between the IRV and the SRV distributions, but 
the former is based on relatively few objects. A compar- 
ison with the C-type IRVs and SRVs, for which the me- 
dian is 9.5 km s^^ (39 objects) and where the fraction of 
low-velocity sources is much lower, indicates that in this 
respect there is a difference between the chemistries. A C- 
type chemistry produces higher gas expansion velocities. 
The large fraction of low- velocity sources in our sample is 



further discussed in Sect 6.5 



6.4. Mass loss and envelope kinematics 

There are two main characteristics of the mass loss pro- 
cess, the stellar mass loss rate and the circumstellar gas 
expansion velocity. The former is to a large extent deter- 
mined by the conditions at the transonic point, i.e., the 
density at the point where the gas velocity goes from being 
sub- to supersonic, while the latter is determined by the 
acceleration beyond this point. Hence, these two proper- 
ties do not necessarily correlate with each other. However, 



in a study of a dust-driven wind Habing et al. (1994) found 



that the relation A/ocuc should apply in the optically thin 



limit. Solving the same set of equations Elitzur & Ivczic 



(2001) found that the dependence becomes even stronger, 
M(xv^, when the effect of gravity is negligible. Therefore, a 
comparison between the two mass loss characteristics may 
provide important results, which any mass loss mechanism 
model must be able to explain. 

In Fig. H we present the mass loss rates and gas expan- 
sion velocities for our sample. There is definitely a trend 



Fig. 10. C0(J=2^1) and C0(J=3^1) spectra (his- 
tograms) obtained at the SEST and the JCMT, respec- 
tively, and model line profiles (thin, solid lines) for the 
low gas expansion velocity source Pup (the beam size 
is given in each panel) 



since the detection rates are high for both samples, we con- 
clude that there is clearly a difference between the mass 
loss properties of AGB-stars of M- and C-type in this re- 
spect. 

Of particular interest for further studies are the 
sources with gas expansion velocities lower than 3kms~^: 
V584 Aql, T Ari, BI Car, RX Lac, and Pup. Such a low 
velocity corresponds to the escape velocity at a distance of 
1.5x10^^ cm for a 1 M© star, i.e., a distance correspond- 
ing to about 100 stellar radii, considerably larger than the 
normally accepted acceleration zone, about 20 stellar radii 
( Habing et al. 1994 ) . In a detailed study of state of the 
art mass loss models Winters et al. (2d00| ) concluded that 
for low radiative acceleration efhciencies only low mass 
loss rates (<3x 10~^ M© yr~^) and low gas expansion ve- 
locities (<5kms~^) are produced (their class B models). 
In these models the gas expands at a relatively constant 
and low velocity beyond a few stellar radii, and it finally 
exceeds the escape velocity at large radii. This can pro- 
vide an explanation for the low velocity sources, but also 
more complicated geometries/kinematics may play an im- 
portant role (Sect. ^.7| ). 

Jura et al. (2002) presented a study of Pup where 
they used the mid-IR cameras on the Keck telescope to 
partly resolve the dust emission. No clear geometrical 
structure is evident, but they derive a dust mass loss rate 
of 10~^ Mq yr~^, which combined with our gas mass loss 
rate, results in a dust-to-gas mass ratio as high as 0.05, 
i.e., about 15 times higher than our average estimate from 
the /i-parameter (Sect [5.21) . This suggests that there are 
problems with the dust and/or the gas mass loss rate es- 
timates of this star. In Fig. ^ we show the model line 
profiles superimposed on our highest quality spectra of 
this object. A very good fit is obtained for a mass loss rate 
of 2.2 X 10"^ M© yr~^, an expansion velocity of 2.1 kms^^, 
and a turbulent velocity width of 1.0 km s^^ (this is higher 
than the 0.5 km s^^ adopted in the modelling of the whole 
sample). However, this object is so extreme that not too 
much faith should be put in the model results despite the 
successful fit. 



Fig. 11. Observed and modelled 
line profiles for R Dor. Upper pan- 
els show the results of the nom- 
inal model of R Dor (_D=46pc, 
M=1.3xlO"''MQyr-\ h=0.7). The 
middle panels show the results for a 
distance three times larger (D=150pc, 
M=5xlO"^MQyr-\ h=1.5). The 
lower panels show the results for a 
photodissociation radius three times 
smaller than that of the nominal model 
(M=3xlO-''M0yr-\ ft=0.05) 



6.6. The line profiles 

The line profiles carry information which was not used to 
constrain the derived mass loss rates. In this section we 
briefly discuss how well the model line profiles reproduce 
the observed ones. In general, the results are satisfactory. 
However, although there are both too double-peaked and 
too rounded model profiles when compared with the ob- 
served ones, there is a trend of too double-peaked (or too 
flat) model line profiles, indeed observed double-peaked 
line profiles are very rare. The discrepancy is worst for 
the J =1^0 line, where double-peaked model line profiles 
are common. There are two possibilities for such a dis- 
crepancy, either the angular size of the emitting region is 
too large in the model or there is an effect due to maser 
emission, which is not reproduced in nature. The former 
can be due to systematically too small distances to the 
sources or too large photodissociation radii. The latter is 
a possibility because for these low mass loss rate objects 



tion are not expected. The double-peakedness is rather 
due to the large angular extent of the emission in the 
model, i.e., the emission is clearly spatially resolved. This 
can be solved by increasing the distance. The somewhat 
uncertain Hipparcos distance of R Dor is 61 pc, but a 
change to this distance leads only to marginal improve- 
ments in the model fitting. We have to increase the dis- 
tance by a factor of three to get acceptable fits to the 
observed data (i?=150pc for which the best-fit results 
are M=5xlO~''M0 yr~^ and h—1.5; note that the de- 
rived temperature structure depends on the distance to 
the source due to the emission being spatially resolved). 
Such a large distance is not obviously compatible with the 
Hipparcos data. Alternatively, we can artificially lower the 
size of the CO envelope b y a factor of three co mpared to 
that given by the model of Mamon et al. (198^ ) [for which 
the best fit results are M=3x 10"'^ Mq yr^^ and /i=0.05; 
for this mass loss rate the photodissociation radius is ac- 



radiatii^e excitation plays an important role and it tends 
to invert preferentially the lower J-transitions. However, 
among the objects with discrepancies there are about as 
many without maser action as with maser action (in the 
model). 

The by far worst discrepancy is obtained for R Dor, 
where all the model line profiles are strongly double- 
peaked. Fig. |ll|. All three transitions are inverted (the 
J=1^0 line over the whole radial range, but the J=2^1 
and J=3— >1 lines only over a part of it), but the optical 
depths are so low that substantial effects of maser ac- 



tually five times larger according to the model of Mamon 
et al. (1988)]. Actual tests of the predictions of the model 



of Mamon et al. (1988) have mainly been done for high 
mass loss rate C-stars ( [Schoier fc Olofsson 2000| , ^00l\ ). 
The model has passed these tests, and it is therefore ques- 
tionable whether it gives results off by a factor of five at 
lower mass loss rates. A possibility exist that the inter- 
stellar UV field is exceptionally strong and hence limits 
severely the CO envelope of R Dor. Alternatively, the mass 
loss is highly variable, and the small outer radius of the 
CO envelope is a consequence of a recent higher mass loss 
epoch. We must conclude that presently the reason for 
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Table 4. Source parameters and model results for those objects with double component line profiles 



Source 


V£Lr. 

type 


P 

[days] 


u 

Inn] 


Ibb 


comp. 


A/f 

[iu iviQ yr J 


1 Kill S 1 


iiu cini 


n 


^2 

A.rcd 


IN 




orlD 


OO 


1 /in 


zzUU 


— r J 

broad 

narrow 




0.3 


O 

1.0 


z.o 
1.1 






1 
1 

1 


RV Boo 


SRb 


137 


280 


2760 


broad 


2.0 


7.0 


1.8 


0.1 


0.7 


4 












narrow 


0.3 


2.3 


0.8 


0.05 


9.9 


4 


X Her 


SRb 


95 


140 


2490 


broad 


1.5 


6.5 


1.5 


0.2 


4.6 


3 












narrow 


0.4 


2.2 


1.0 


0.03 


4.5 


3 


SV Psc 


SRb 


102 


380 


2450 


broad 


3 


9.5 


1.9 


0.1 


2.1 


3 












narrow 


0.3 


1.5 


1.0 


0.05 


20.7 


3 



the major discrepancy between the observational data of 
R Dor and our modelling results is not clear. A complicat- 
ing factor is that the J— 1-^0 line profile is time variable 
(Nyman, priv. comm.), and that it at times looks rather 



peculiar [compare e.g. the spectrum shown in Lindqvist 
et al. (|1992D]. 



Finally, we note here that Olofsson et al. (1993) in their 



sample of C-type SRVs and IRVs found five sources (about 
10%) with remarkable CO radio line profiles, sharply 
double-peaked with only very weak emission at the sys- 
temic velocity. These were s ubsequently shown t o origi- 
nate in large, detached shell s (Olofsson et al. 1996), which 
are ge omet rically very thin ( Lindqvist et al. 1999 ; [Olofsson 
et al. 2 000 ), presumably an effect of highly episodic mass 
loss. We have found no such source in our sample of M- 
stars. Hence, in this respect there must be a difference in 
the mass loss properties between the two chemistries. 

6. 7. Double-component line profiles 

Knapp et al. (1998) and [Kerschbaum fc Olofsson (1999|) 
were the first to discuss more thoroughly the small num- 
ber of objects with line profiles which can be clearly di- 
vided into two components, a narrow feature centred on 
a broader pla teau. Such sources exist both among M- and 
C-type stars ( Knapp et al. 1998). The origi n of such a line 
profile is still not clear. Knapp et al. (19981) argued that it 
is an effect of episodic mass loss with highly varying gas 
expansion velocities. Alternatively, it can be an effect of 
complicated geometries/ki nematics. The first spa tial in- 
formation was provided by Kahane fc Jura (1996 ). A CO 



radio line map towards the M-star X Her suggested that 
the broad plateau is a bipolar outflow, while the narrow 
feature was not spatially resolved. Bergman et al. (2000 ) 
produced CO radio line interferometry maps of the M-star 
RV Boo. In this case the brightness distributions suggest 
that the broad plateau emission comes from a circumstel- 



lar disk with Keplerian rotation. Kahane et al. (1998) and 
Jura fc Kahane (199S| ) interpret the narrow CO radio line 
features which they observe in a few cases as originating 
in reservoirs of orbiting gas (these sources do not have 
distinct double-component line profiles). It is fair to say 
that no consensus has been reached about these peculiar 
circumstellar emissions. 



In our sample we have four sources of this type, 
EP Aqr, RV Boo, X Her, and SV Psc, all of them SRVs. 
We have simply decomposed the emission into two com- 
ponents for each source, assuming that the emissions are 
additive. Mass loss rates and gas expansion velocities were 
determined in the same way as for the rest of our objects. 
This is probably a highly questionable approach for both 
components. The results, as well as some source informa- 
tion, are given in Table ^. Knapp et al. (19"98[ ) derived 
mass loss rates for EP Aqr and X Her which are within 
a factor of two of our estimates (both for the narrow and 
the broad components). 

Not unexpectedly the mass loss rates are higher for 
the broader component by, on average, an order of magni- 
tude. The fits to the narrow components result in very low 
gas expansion velocities. Indeed, so low, e.g., l.Okms"^ in 
the case of EP Aqr, that an interpretation in the form of 
a spherical outflow is put to question. On the other hand, 
the relation between mass loss rate and gas expansion ve- 
locity is the same for these objects as for the rest of the 
sources. Fig. ^ The narrow component gas appears cooler 
than the broad emission gas in those three cases where an 
h can be estimated. This may be an accidental result, but 
it can also provide a clue to the interpretation. 



7. Conclusions 

We have determined mass loss rates and gas expansion 
velocities for a sample of 69 M-type IRVs (22 objects) and 
SRVs (47 objects) on the AGB using a radiative trans- 
fer code to model their circumstellar CO radio line emis- 
sion. We believe that this sample is representative for 
the mass losing stars of this type. The uncertainties in 
the estimated mass loss rates are rather low within the 
adopted stellar/circumstellar model, typically less than 
±50%. However, a sensitivity analysis shows that for these 
low mass loss rate stars there is a considerably uncertainty 
due to the stellar luminosity, the size of the CO enve- 
lope, the CO abundance, and as usual the distance to the 
source. We find that the mass loss rates determined by the 
detailed radiative transfer analysis differ by almost an or- 
der of magnitude from those obtained by published mass 
loss rate formulae. 
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The (molecular hydrogen) mass loss rate distribution 
has a median value of 2.OxlO~''M0 yr~^, and a minimum 
of 2.0xl0~* Mq yr~^ and a maximum of 8x10"'^ yr~-^. 
M-type IRVs and SRVs with a mass loss rate in excess of 
5x10~''Mq yr~^ must be very rare, and in this respect 
the regularity and amplitude of the pulsation plays an 
important role. We also find that among these mass losing 
stars the number of sources with mass loss rates below a 
few 10~* M0 yr"-'^ must be small. 

We find no significant difi'erence between the IRVs 
and the SRVs in terms of their mass loss characteristics. 
Among the SRVs the mass loss rate shows no dependence 
on the period. Thus, for these non-regular, low-amplitude 
pulsators it appears that the pulsational pattern plays no 
role for the mass loss efficiency. 

We have determined temperatures for our sample stars 
by fitting blackbody curves to their spectral energy distri- 
butions. These blackbody temperatures have been shown 
to correlate reasonably well with the stellar efi'ective tem- 
peratures. The mass loss rates of our stars show no corre- 
lation at all with these stellar blackbody temperatures. 

The gas expansion velocity distribution has a median 
of 7.0kms~^, and a minimum of 2.2kms~^ and a maxi- 
mum of 14.4kms~^. No doubt, these objects sample the 
low gas expansion velocity end of AGB winds. The fraction 
of objects with low gas expansion velocities is high, about 
30% have velocities lower than 5kms~^. There are four 
objects with gas expansion velocities lower than 3kms^^: 
V584 Aql, T Ari, BI Car, RX Lac, and Pup. These 
objects certainly deserve further study. 

We find that the mass loss rate and the gas expansion 
velocity correlate well, M(xvl^^, even though for a given 
velocity (which is well determined) the mass loss rate may 
take on a value within a range of a factor of five (the 
uncertainty in the mass loss rate estimate is lower than 
this within the adopted circumstellar model). The result 
is in line with theoretical predictions for an optically thin, 
dust-driven wind. 

A more detailed test of the CO modelling is provided 
by the shape of the line profiles. In general, the fits are 
acceptable, but there is a trend that the model profiles, in 
particular the J=1^0 ones, are more flat-topped, or even 
weakly double-peaked, than the observed ones. An excep- 
tional case is R Dor, where the high-quality, observed line 
profiles are essentially flat-topped, while the model ones 
are sharply double-peaked. Acceptable fits are obtained 
by increasing the distance to the star or by artificially 
decreasing the size of the CO envelope. 

The sample contains four sources with distinctly 
double-component CO line profiles: EP Aqr, RV Boo, 
X Her, and SV Psc (all SRVs). We have modelled the 
two components separately for each star and derive mass 
loss rates and gas expansion velocities using the same cir- 
cumstellar model as for the rest of the sample. The result- 
ing mass loss rates and gas expansion velocities show the 
same positive correlation as that of the other objects. At 
present, the exact nature(s) of these objects is unknown. 



We have compared the results of this M-star sample 
with a similar C-star sample. The mass loss rate distribu- 
tions are comparable, suggesting no dependence on chem- 
istry for these types of objects. Likewise, the mass loss 
rates of the C-stars show no correlation with stellar tem- 
perature or period. The gas expansion velocity distribu- 
tions though are clearly different. The fraction of low ve- 
locity sources is much higher in the M-star sample. In both 
cases there is a correlation between mass loss rate and gas 
expansion velocity, although the detailed relations are dif- 
ferent. Our crude estimates of the dust properties, through 
the gas-grain collision heating term, indicate that the two 
samples have similar gas-to-dust ratios and that these dif- 
fer significantly from that of high mass loss rate C-stars. 
This also means that the gas-CSEs due to low mass loss 
rates are cooler than expected from a simple extrapolation 
of the results for IRC-l-10216. Finally, we find no example 
of the sharply double-peaked CO line profile, which is ev- 
idence of a large, detached CO-shell, among the M-stars. 
About 10% of the C-stars show this phenomenon. 
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